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The study of quantum computation has been motivated by the hope of finding efficient quantum 
algorithms for solving classically hard problems. In this context, quantum algorithms by local 
adiabatic evolution have been shown to solve an unstructured search problem with a quadratic 
speed-up over a classical search, just as Grover's algorithm. In this paper, we study how the 
structure of the search problem may be exploited to further improve the efficiency of these quantum 
adiabatic algorithms. We show that by nesting a partial search over a reduced set of variables into a 
global search, it is possible to devise quantum adiabatic algorithms with a complexity that, although 
still exponential, grows with a reduced order in the problem size. 



I. INTRODUCTION 

Grover's quantum algorithm solves an unstructured 
search problem in a time of order V^V, where N is the 
dimension of the search space, which corresponds to a 
quadratic speed-up over a classical search 0. This algo- 
rithm is proved to be optimal in the case of unstructured 
search problems Q- Naturally, it can also be used to 
solve a structured search problem with a quadratic speed- 
up over a naive classical search that would exhaustively 
check every possible solution. However, exploiting the 
structure of the problem is well known to lead to bet- 
ter classical search algorithms. It is therefore tempting 
to imagine that better quantum search algorithms may 
be devised similarly by exploiting the problem structure. 
Following this, Cerf, Grover and Williams have shown 
that this could be done by partitioning the unknown 
variables into two (or more) sets and nesting a quan- 
tum search over one set into another search over two (or 
more) sets, yielding an average complexity of order ^/N a , 
with a<l| 

While this algorithm, as well as Grover's original al- 
gorithm, stay within the standard paradigm of quantum 
computation based on quantum circuits, a new type of 
quantum algorithms based on adiabatic evolution has 
been introduced lately 0, Q. In particular, a quantum 
adiabatic analogue of Grover's search algorithm has been 
proposed in Ref. 6] , which works for unstructured search 
problems. The use of quantum adiabatic algorithms has 
also been analyzed for solving structured problems such 
as k-SAT, but in such a way that until now only a numer- 
ical study has been possible Q. Recently, the study of 
quantum adiabatic algorithms has progressed even fur- 
ther after Aharonov and Ta-Shma demonstrated their 
universality, namely that any problem that may be solved 
efficiently in the standard (circuit-based) model of quan- 
tum computation can also be solved efficiently by a quan- 
tum adiabatic algorithm 8]. This provides a strong in- 
centive to the search for new quantum adiabatic algo- 
rithms. 

The purpose of this paper is to bring the ideas of nested 
quantum search and quantum adiabatic computation to- 



gether, in order to devise a new quantum adiabatic algo- 
rithm adapted to structured problems. More specifically, 
we will show that an adiabatic search over a subset of the 
variables can be used to build a better initial Hamilto- 
nian for the global adiabatic search. With this adiabatic 
algorithm, we recover the same complexity as the nested 
circuit-based algorithm of Ref. Q . 



II. ADIABATIC THEOREM 

Let us briefly recall the Adiabatic Theorem and how it 
may be used to design quantum algorithms by adiabatic 
evolution. 

We know that if a quantum system is prepared in the 
ground state of a time-independent driving Hamiltonian, 
it remains in this state. The Adiabatic Theorem states 
that if this Hamiltonian becomes time-dependent, the 
system will still stay close to its instantaneous ground 
state as long as the variation is slow enough. 

More specifically, if \E ;t) and \Ei;t) are the ground 
and first excited states of the Hamiltonian H(t), with 
energies Eo(t) and E\(t), we define the minimum gap 
between these eigenvalues 



min [E^t) - E (t)} 

0<t<T 



(1) 



and the maximum value of the matrix element of dH/dt 
between the eigenstates as 



max 

0<*<T 



,dH_, 

1 dt ' 



(2) 



with (^r)i,o = (-Ei;i|^jr|.Eo;*)' The Adiabatic Theorem 
states that if we prepare the system at time t = in its 
ground state \Eq;0) and let it evolve under the Hamilto- 
nian H(t), then 
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where e<l. 

Now, we may apply to the system a Hamiltonian for 
which the ground state encodes the unknown solution 
of a problem. According to the Adiabatic Theorem, we 
know that we may get the system very close to this solu- 
tion state by preparing it in the (known) ground state of 
another Hamiltonian, and then by progressively chang- 
ing it to the Hamiltonian of our problem. This simple 
idea is central to the quantum algorithms by adiabatic 
evolution 0,0. 



III. QUANTUM SEARCH BY LOCAL 
ADIABATIC EVOLUTION 



As exposed in [£g , this principle may be used to perform 
a quantum search. Suppose that among N states, we 
have to find the M-times degenerate ground state of a 
Hamiltonian 



H f = i- WW> 



(5) 



where M. is the ensemble of solutions (of size M). We 
initially prepare the system in an equal superposition of 
all could-be solutions: 



(6) 



This superposition is the ground state of the following 
Hamiltonian: 



Hi=I-\s)(s\. 



(7) 



We now apply Hi to the system and switch adiabatically 
to Hf. If we perform an adiabatic evolution 



H(t) = (l-8(t))Hi + 8(t)H f , 



(8) 



where s(t) is a (carefully chosen) monotonic function with 
s(0) = and s(T) = 1, wc will finally obtain a state close 
to a ground state of Hf. 



lV>/> 
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as long as 



T = 




(9) 



(10) 



This algorithm is refered to as local because s(t) is chosen 
such that the Adiabatic Theorem is obeyed locally, at 
each time (see [g for details). 



each additional variable considered, new constraints ap- 
pear and reduce the set of satisfying assignments. This 
corresponds to most problems encountered in practice (k- 
SAT, graph coloring, planning, combinatorial optimiza- 
tion, . . . ). 

For a set of ha variables denoted as A, there is a corre- 
sponding set of constraints Ca ■ We may define a function 
Ja that tells if an assignment of the variables in A satis- 
fies the constraints in Ca- 



f A ■ (z d ) n * -> {0,1} 

if a; does not satisfy Ca 

1 if x satisfies Ca, 



(11) 



where d is the number of possible assignments for each 
variable (d = 2 for bits). As quantum gates have to be 
reversible, the quantum equivalent of this function will 
be an oracle: 



o A ■n NA ®n 2 ^'H 



N A < 



m 2 



\x)®\y®fA(x)), 

(12) 

where Na = d nA . It is shown in Ref. |2( that this oracle 
is closely related to a Hamiltonian whose ground states, 
of energy 0, are the basis states encoding a satisfying 
assignment and whose excited states, of energy 1, are all 
other basis states: 



H A \x) 



\x) li x does not satisfy Ca 
if x satisfies Ca 



or 



H a = Ia- E \ m A){mA\ 

ttia^Ma 



(13) 



(14) 



where Ma is the set of satisfying assignments for the 
variables in A. It is possible to reproduce the application 
of this Hamiltonian during a certain time by using two 
calls to the quantum oracle Oa (see @ for details). 

Now suppose we consider a larger set of variables 
71 ab = fiA + n-B that have to satisfy a set of constraints 
Cab D Ca- To discriminate between assignments satis- 
fying Cab or not, we will use an oracle Oab or a cor- 
responding Hamiltonian Hab defined as in Eqs. Ijl2|l 
and i|13|) . The basic idea of our structured search will be 
to find the solutions of Cab by first building the assign- 
ments of the riA primary variables satisfying Ca, then by 
completing them with all possible assignments of the tib 
secondary variables, and finally by searching among these 
could-be solutions the global satisfying assignments. 



STRUCTURED SEARCH WITH TWO 
LEVELS OF NESTING 



IV. STRUCTURED PROBLEMS 

In this article, we consider a class of problems where 
one has to find an assignment for a set of variables. For 



This problem is of the same type as the one consid- 
ered in |3J, for which the technique of nesting was intro- 
duced in the context of the traditional implementation 
of Grover's algorithm on a conventional quantum circuit. 
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Here, we apply this technique to the adiabatic quantum 
search algorithm. 

Suppose we divide the variables of our problem into 
two subsets A (ua elements) and B (ns elements). First, 
we will perform a search on the variables in A using the 
Hamiltonian Ha that encodes the constraints in C A : 



H a =Ia- E \m A )(m A \ 

m A GMA 



(15) 



Then we will use the Hamiltonian Hab acting on all vari- 
ables in A U B and encoding the whole set of constraints 
Cab 



Hab = Iai 



(16) 



E 



\mA){m A \ <8> \m B ){rn B \ 



(m A ,m B )£M A i 



to construct a superposition of the solutions of the full 
problem Mab- A final measurement of the quantum 
register then gives one of the global solutions at random. 



A. Adiabatic search on the primary variables 

The preliminary search on the variables in A is a sim- 
ple unstructured search as explained in section 11111 As 
there are ua variables in A, the corresponding Hilbert 
space is of dimension Na = d nA . Let Ma be the number 
of solutions in Ma- Performing an adiabatic quantum 
search, we may thus transform the initial state 



1 



E 



(17) 



into a state close to the uniform superposition of all so- 
lutions in Ma 



in a time of order 




(18) 



(19) 



B. Adiabatic search on the secondary variables 

We will now perform a preliminary search in the 
Hilbert space of dimension Nb = d riB of the secondary 
variables in B by extending the partial solutions Im^)- 
We prepare the variables in B in a state that is the uni- 
form superposition 



SB 



3)B. 



(20) 



Globally, the system is thus in the superposition: 



AB 



\lpr, 



SB 



1 



J2 \ m *) ® \j)b, (21) 



-> A eM A 



where some terms correspond to a global solution of the 
problem [{m A ,j) € Mab satisfying all constraints in 
Cab) and the others to a partial solution only [rriA € -Ma 
satisfies Ca but (rriA,j) £ Mab does not satisfy Cab]- 
We now divide the set Ad a of solutions of Ca into two 
subsets: M A will be the set of m^'s for which there exists 
at least one solution (m A , ws) of Cab and M A the set 
of m^'s for which there is no such solution. 

M S A = {m A € M A \3m B ,(m A ,m B ) e M AB }{22) 
Mf = {m^^lVj.Kj)^^} (23) 



M S A U M A S . We may 
now rewrite our initial state (|21|l as 



Of course, we thus have M A 

j 
i 



\4>o)ab = 



VMaNe 



1 



VM A N B 



\ m A) ® \])b 

M A S 

E \ m A)®\j)B- 



j SAT 



(24) 



a A £M-; 



In the first part of this expression, no term correspond 
to a solution of the full problem, whereas in the second 
part, some terms do and others do not. The goal of this 
stage of the computation will be to increase the ampli- 
tude of the solution terms in this last part. To achieve 
this, we perform an adiabatic evolution using as initial 
Hamiltonian 



Hi=I A ® (I B ~ \s b )(sb\), 



(25) 



that has \ipo) ab as a ground state. The final Hamiltonian 
will be 



H f = Hab - H a ®Ib- 



(26) 



We see that these Hamiltonians share the following prop- 
erties: 

1. They do not act on states \i)a®\sb) corresponding 
to assignments i of Af A that do not satisfy C A : 
Hi,f\i)A®\s B ) =0 Vi<£M A - 

2. They do not couple states corresponding to differ- 
ent rriAS: 

B (j\ ® (m A \Hij\m' A ) <g> \j') B = V m A ^ m! A e 
M A , Vj,j'£JV b . 

It follows that the instantaneous Hamiltonian of the adia- 
batic evolution H(t) satisfies the same properties. Keep- 
ing this in mind, it may easily be shown that the effect 
of the adiabatic evolution will be to perform independent 
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adiabatic searches for each tua € Ma- More precisely, 
each term in |Vo)ab 



7= E l TOA > \fi B 

'iVR 777 



(27) 



will evolve to 



mA m B eM B/r , 



\m A ) <8> |m B ), (28) 



as long as 



M 



a B/m A 



(29) 



where M.B/m A i s the set of ma's such that (mx, ms) S 
jVtyiB and M B / mA is the number of these elements. For 
this condition to be satisfied for all m^'s simultaneously, 
we must take 



T B = max T„ lA = O 

m A 



N B 



min mA M B / 



(30) 



At the end of this second stage, we thus have constructed 
a state close to 



AB) 



1 V e^A\ mA ) 



m A eM% 



E 



\m B ) 



mA m B GM B/r , 



where 4> mA 's are phases appearing during the evolution, 
1 



E I"»a)®U)b (32) 



M A N B mA€M »s 

ieU B 

1 £ e ^| TOA ) 



E 



V M B/m A m B €-M B/mjl 



fib; 



(33) 



and (resp. Mf) is the number of elements in set 

M^ s (resp. M^). 



C. Global adiabatic search 

Stages A and B define a unitary evolution U that ap- 
plies the initial state \sa) ® \sb) onto IVab): 



C/|sa) ® \s B ) 



>AB) 



(34) 




(35) 



In this state, we now need to decrease the amplitude of 
the first term, corresponding to partial solutions only, 
and increase the amplitude of the second term, corre- 
sponding to global solutions. This could be realized effi- 
ciently by performing an adiabatic search using as initial 
Hamiltonian: 

Hi = Iab-\iPab){^ab\ (36) 
« U{Iab~\sa){sa\®\sb){sb\)U ] (37) 
k, UH Q U\ (38) 

where H = Iab~ \sa){sa\®\sb){sb\j an d as final Hamil- 
tonian 



H 



f 



H AB 
Iab - 



(39) 



\mA)(m A \ ® \m B )(mB\ 

(m A ,m B )GM AB 



during a time 



Tn = 



! M A 
Ml 



(40) 



Unfortunately, we do not have access to Hi, so that the 
interpolating Hamiltonian H(s) = (l — s)Hi+sHf cannot 
be applied directly. However, we will see in Section IVT1 
that the basic steps of the quantum circuit implementa- 
tion of this adiabatic algorithm only require the applica- 
tion of Hi during a particular time i, that is 



-iHit 



-iUHoUh _ 



Ue~ lHot U^. 



(41) 



Hence, each application of Hi during a time t will 
be equivalent to sequentially applying U>, e~ lH<)t , and 
U, which means performing the adiabatic evolution U 
(stages A and B) backwards, then applying H for a time 
t, and finally rerun U forwards (stages A and B). 

In section fvTI we will see that, when discretizing the 
evolution, we must take a number of steps rc of order Tc- 
We may now evaluate the complexity of this algorithm. 
As it consists of tq steps, each involving two applica- 
tions of U or f/T , that last a time of order Ta + T B , the 
algorithm finally takes a time of order: 




m A M B l JrlA 
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Let us notice that with the same hypothesis as in Ref. Q , 
namely 



M 



B/m A 



1 V m A , 



then Mf = Mabi an d the computation time is 



O 



AB 



(44) 



(45) 



i.e., the same complexity as the equivalent circuit-based 
algorithm exposed in Ref. A more detailed analysis 
of this complexity will be performed in section IVIII 



VI. DISCRETIZING THE ADIABATIC 
EVOLUTION 

A. General method 

The method to implement a global adiabatic evolution 
algorithm on a discrete quantum circuit was initially ex- 
posed in [5j , and was extended to the case of a local adi- 
abatic evolution algorithm in |j. Let us quickly review 
this method, that uses two successive approximations. 

The first approximation consists in cutting the evo- 
lution time T in r intervals AT = — and replacing the 
continuously varying Hamiltonian H (t) by a Hamiltonian 
H'(t) that is constant during each interval AT and varies 
at times tj = jAT only. 



H'(t) = H(tj) if tj-i <t<tj 



(46) 



It is shown in |9j that for H(s) = (1 — s)Hi + sHf with 
s = s(t), this approximation introduces a global error for 
the evolution 



lUW-U'WfoK^jlHi-Hfh, (47) 

where ||A|||2 = rcLaxii i^v ||=i ll^l^)!! is the operator norm of 
A. Our algorithm now requires r steps of the form 

where Sj = s(tj). As we are able to apply H\ and Ht sep- 
arately but not necessarily a simultaneous combination 
of them, we will approximate C/j by 



(J" = e -t(l-« J )H 4 AT e -i* J H / Ar < 



This will result in an error 



\\u'(T)-]Ju;%eO(^Wi,.ii nil2 , 



(49) 



(50) 



B. Application to a structured quantum search 

We now consider the case of a structured quantum 
search. We could apply the discretization procedure to 
all three stages (A, B, C) of our algorithm in order to im- 
plement it on a quantum circuit, but we will concentrate 
on stage C which is the only one that requires discretiza- 
tion. Nonetheless, it is easy to show that stage A (resp. 
B) would require a number of steps ta (resp. re) of the 
same order as the computation time Ta (resp. Tg). 

For the final stage, the global adiabatic search, the 
Hamiltonians Hi and Hf are defined in Eq. (|36(l - (|39|l . 
Evaluating the errors introduced by the approximations, 
we find: 



\\H, - H 



/III 2 



< 1 



Wi,H f ]\h < 1 
and, asT c = 0(y^), 

l\U(T)-U'(T)\j 2 G O 

\tu'(T)- n^'Ba e O 



'Ml 

M A 



(51) 
(52) 



\ 



(53) 



(54) 



Therefore, as announced in section^] we have to cut our 
evolution in a number of steps rc — O (^Jm&J °f the 
same order as Tq- Each step j will take the form: 



jjii = e -i(l-s 3 )H^T e -i Sj H f AT 

w Ue~^ 1 ' Sj)HoAT U t e" iSjHf 



AT 



(55) 
(56) 



(see [3 for details). 



where the applications of Hamiltonians Hq during a time 
(1 — Sj)AT and Hf during a time SjAT may be realized 
by the procedure exposed in 0. 



VII. COMPLEXITY ANALYSIS 

To estimate the efficiency of this algorithm, we will fol- 
low the same development as in |3J: as we have seen in 
Section |VJ under the assumption l|44l) that we will con- 
sider here, the complexity of this adiabatic algorithm has 
exactly the same form as its circuit-based counterpart. 

First of all let us define a few concepts (for details 
here and throughout this section, we refer the reader to 
Ref. 0). The structured search problem is to find an 
assignment of hab = n-A + nB variables among d pos- 
sibilities and satisfying e constraints, each involving at 
most k of these variables. We define as a ground instance 
an assignment of all the variables involved in a particular 
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constraint. A ground instance will be said to be nogood 
if it violates the constraint. Let £ be the number of those 
no-good ground instances. 

Empirical studies show that the difficulty of solving a 
structured problem essentially depends on four parame- 
ters: the number of variables uab, the number of pos- 
sible assignment per variable d, the number of variables 
per constraint k, and the number of no-good ground in- 
stances £. Intuitively, we understand that if £ is small, 
there are many assignments satisfying the constraints so 
the problem is easy to solve. On the contrary, if £ is 
large, the problem is over-constrained and it is easy to 
show that there is no solution. More precisely, it may 
be shown that for fixed uab and d, the average diffi- 
culty may be evaluated by the parameter = £,/uab- 
The hard problems will be concentrated around a criti- 
cal value (3 C . 

Let us now estimate the complexity 1)45 p. Let p(n) be 
the probability that a randomly generated assignment of 
the n first variables satisfies all the constraints involv- 
ing these variables. We then have Ma = p(riA)d nA and 
Mab = p{nAB)d nAB while it is shown in [3| that: 



p(n) W d- n AB(0/Pa){n/n AB ) k 

Eq. I|45|) becomes 



T = 



I (Jp-A _|_ ^ ( ]riAB[l-(P/l3< : )(nA/riAB) k \ 



or, with a = V d nAB and x — ua/tiab'- 



T = 



(57) 



(58) 



(59) 



We now optimize x, the fraction of variables for which we 
perform a partial search, to minimize the computation 
time. We have to solve the equation 



0c 



kx 



k-l 



x'+x— 1 



(60) 



which, for large a (that is large uab) approximately re- 
duces to 



k 



1 = 0. 



(61) 



The solution of this equation a (0 < a < 1) corresponds 
to the optimal partial search we may perform such that 
the complexity grows with the smallest power in d for 
n ab ^ oo. This optimal computation time may then be 
written as 



T = 



2a a 




(62) 



Let us now consider the hardest problems for which (3 : 
[3 C . For these problems, the complexity reads: 



T = 



(63) 



that we may immediately compare to the complexity of 
an unstructured quantum search O The gain 

in the exponent a depends on k through Eq. (|61|l For 
instance, we find a = 0.62 for k = 2, a = 0.68 for k = 3 
and a-* 1 when k — > oo. 

To compare these results with a classical algorithm, 
let us consider a specific problem, the satisfiability of 
boolean formulas in conjunctive normal form, or fc-SAT. 
For 3-SAT, which is known to be NP-complete, the best 
known classical algorithm has a worst-case running time 
that scales as O (2 0A5nAB ) 10], while, as a = 0.68 for 
k = 3, our quantum adiabatic algorithm has a compu- 
tation time of order O (2°- 34 ™^ B ), which is a slight im- 
provement. Let us also notice that this scaling could 
be further improved by using several levels of nesting, 
i.e. by replacing the preliminary search over the primary 
variables by a another nested structured search (see the 
analysis of the circuit-based counterpart of this idea in 
the Appendix of |3j). 



VIII. CONCLUSION 

We have introduced a new quantum search algorithm 
combining the approach based on local adiabatic evo- 
lution developed in |6J and the nesting technique intro- 
duced in It allows one to adiabatically solve struc- 
tured search problems with an improved complexity over 
a naive adiabatic search that would not exploit the struc- 
ture of the problem. 

The basic idea is to perform a preliminary adiabatic 
search over a reduced number of variables of the problem 
in order to keep only a superposition of the assignments 
that respect the constraints of this partial problem, and 
then to complete these partial solutions by finding satis- 
fying assignments for the remaining variables. We have 
seen that to implement this algorithm, the global adi- 
abatic evolution (stage C) has to be discretized, which 
makes it possible to nest the preliminary search (stages 
A and B) into the global one. Each step of the algorithm 
requires to alternate partial adiabatic searches backwards 
and forwards with global search operations. 

A complexity analysis shows that the average compu- 
tation time of this adiabatic algorithm, although still 
exponential, grows with a reduced exponent compared 
to quantum unstructured search algorithms or classical 
known algorithms to solve a problem such as 3-SAT. 
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